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We introduce, and numerically study, a system of N symplectically and globally coupled standard 
maps localized in a d = 1 lattice array. The global coupling is modulated through a factor r~ a , being 
r the distance between maps. Thus, interactions are long-range (nonintegrable) when < a < 1, 
and short-range (integrable) when a > 1. We verify that the largest Lyapunov exponent Am scales 
as Am oc N~ K( - a \ where n(a) is positive when interactions are long-range, yielding weak chaos 
in the thermodynamic limit N — » oo (hence Am -» 0). In the short-range case, k(cj) appears to 
vanish, and the behaviour corresponds to strong chaos. We show that, for certain values of the 
control parameters of the system, long-lasting metastable states can be present. Their duration t c 
scales as t c oc N 13 ^, where /3(a) appears to be numerically consistent with the following behavior: 
f3 > for < a < 1, and zero for a > 1. All these results exhibit major conjectures formulated 
within nonextensive statistical mechanics (NSM). Moreover, they exhibit strong similarity between 
the present discrete-time system, and the a-XY Hamiltonian ferromagnetic model, also studied in 
the frame of NSM. 
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I. INTRODUCTION 

The study of complex systems has been lately one the 
most active areas of investigation. This multidisciplinary 
subject has interest in natural sciences as well as in social 
and artificial systems. In physics, a major area involved 
is statistical mechanics. It is precisely for such systems 
that nonextensive statistical mechanics (NSM) [l| ap- 
pears to have its applicability. Typical features present in 
complex systems are long-range interactions, long-term 
memory, fractal phase-space structure, scale-free network 
structure, or even combinations of these characteristics. 
Quite frequently these systems cannot be correctly de- 
scribed by the well-established Boltzmann-Gibbs statis- 
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tical mechanics (BGSM). Indeed, they often fail to ver- 
ify its basic assumptions (equiprobability of phase-space 
occupation Q and ergodicity Q at the thermal equilib- 
rium state). NSM generalises the usual BGSM formalism 
through the definition of the q-entropy S q = jgMg)] 

(with Si = Sbg = — J dxp(x) lnp(x), where BG stands 
for Boltzmann-Gibbs) under simple constraints Q. 

This theory has been extensively applied to nonlinear 
dynamical systems 0, IS IE 0> as they are one of the 
most widespread and useful ways for modelling complex 
phenomena. For example, Hamiltonian systems are at 
the core of physics and consequently their relevance is 
evident. There has been an important number of results 
indicating that certain conservative models (e.g., the a- 
XY model H 0, the a-Heisenberg model [lfllljj, the 
a-Lennard- Jones-like gas |l2j) can present a behaviour 
departing from the one predicted by the BG formalism. 
For certain classes of initial conditions and parameters, 
the system is prevented from virtually attaining its ex- 
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pected equilibrium state in finite time when N — > oo. 
In other words, if the N — ► oo limit is taken before 
the i — ► oo limit, the system becomes nonergodic. This 
nonergodicity is also reflected in a variety of anomalous 
behaviours such as non-Gaussian momenta probability 
density functions, negative specific heat, aging, and oth- 
ers 0, 0, 0, 0. On the other hand, other (sim- 
pler) nonlinear dynamical systems, as for example maps, 
have emerged in many contexts, very often exhibiting 
new and interesting results. Among these systems we 
find the logistic map, the standard map, coupled map 
lattices, discretised Lorentz gas, and many others jl7j . 
Both low- and high-dimensional discrete-time dynamical 
systems have been studied under the framework of NSM, 
exhibit ing various connections with Hamiltonian dynam- 
ics mm HJ. 

It is known quite well that the number of degrees 
of freedom of a system defines the possibilities that its 
dynamics may approach. A clear illustration of this 
fact in low-dimensional systems (e.g., the 2-dimensional 
Chirikov- Taylor standard map) is the KAM tori, a com- 
plex structure in the map phase-space that separates 
chaotic from regular regions. While much is known for 
systems with few degrees of freedom, the situation is 
more intricate when many elements are involved. This is 
because the T phase-space grows in an extremely rapid 
manner. In addition to this, studying continuous-time 
many-body systems can turn out to be very difficult, if 
not impossible at all, due to the considerable computa- 
tional time needed to integrate the evolution equations. 
This is even more so when long-range forces are involved, 
as it is not justified to neglect any interaction between 
the elements. 

An alternative for this problem is the substitution of 
the continuous-time system for simpler discrete-time sys- 
tems, such as maps, which conserve many of the impor- 
tant features of the physics involved. This substitution 
is in fact justified interpreting a dimensional map sys- 
tem as the intersection (in the Gibbs T phase-space) 
of a Poincare plane with the orbit of a higher dimen- 
sional Hamiltonian. More precisely, if we consider a 
time-independent Hamiltonian system with (N + 1) de- 
grees of freedom ( (2 N + 2)-dimensional T phase-space), a 
(2A r )-dimensional symplectic map is the result of taking 
a Poincare section over the constant energy hypersur- 
face |Tt| . The recurrence time is discrete and the map 
is useful in displaying various properties of the original 
Hamiltonian system. 

NSM concepts have been shown to have an important 
role in low-dimensional maps at the edge of chaos. The 
archetypical logistic map (a D = 1, dissipative, nonlin- 
ear map) is one of the most widely studied dissipative 
systems. In part due to its simplicity, it is often used to 
illustrate many of the most important features of chaos. 
In recent works, Robledo and Baldovin [2l| have ana- 
lytically proved, using standard renormalisation-group 
techniques, that the dynamics of the logistic map at 
its critical point is unmistakenly well-described within 



a NSM frame. The sensitivity to initial conditions is 
a g-exponential function |22|. and is related to the en- 
tropy production through the ^-generalised Pesin-like 
identity, linking the sensitivity to initial conditions to 
the g-entropy S q with q = 0.2445 . . .. Moreover, the 
logistic-map with noise (a Langevin-like generalisation 
of the usual logistic map) has been found to present 
two-step relaxation processes and aging (presenting in- 
teresting common points with slow glassy dynamics |23|). 
Other aspects of the NSM formalism were also studied in 
the more general case of one-dimensional z-logistic family 
of maps considering the attractor fractal nature |24| . Re- 
garding two-dimensional (D — 2) dissipative systems (as 
the Henon map or its linearised version, the Lozi map), 
results indicate that it presents the same value of q as 
the logistic map, therefore suggesting a common univer- 
sality class |25j. Regarding two-dimensional (D = 2) con- 
servative maps, a very interesting example, the triangle 
Casati-Prosen map (mixing, ergodic, but with vanishing 
Lyapunov exponent), has been recently been studied in 
connection to the entropy S q with q = f26| . 

Moving further into conservative discrete-time sys- 
tems, two-step relaxation has been also observed at the 
edge of chaos for the Chirikov- Taylor standard map, 
a paradigmatic one for 2-dimensional symplectic maps. 
This map has been thoroughly studied and explained by 
means of the KAM-theorem 17]. This map is known to 
be area preserving, hence its variables are often consid- 
ered as "coordinate" and "momentum" , in analogy with 
Hamiltonian systems. Along this lines, it is appealing to 
look to the more general case of a symplectic system of 
N coupled standard maps (D — 2N,N > 1). Different 
efforts were done for values of N as low as N = 2 up to 

N = 500 mill EH- In these cases > the same two-step 
relaxation was found as well as other features pointing to 
a NSM applicability. But in all these cases, the coupling 
was done trough the momentum variables. Even though 
this has its own interest, it would be instructive to see the 
effects of coupling such a system through the coordinate 
variables, as in more realistic situations. 

In this work we will study a high-dimensional glob- 
ally coupled conservative map system that, as discussed 
above, presents many of the characteristics of Hamilto- 
nian dynamics. Our purpose is to contribute to the un- 
derstanding of the role that NSM plays in the anomalous 
features present in long-range-interacting dynamical sys- 
tems. Our results show interesting similarities between 
this map model and many-body long-range-interacting 
Hamiltonian dynamics, in particular the a-XY model. In 
both cases, long-lasting metastable states, as well as weak 
chaos, are present in the thermodynamic limit. Both fea- 
tures point out that under certain initial conditions, or- 
bits do not visit with equal probability the entire T phase 
space, or in other words there is a failure in ergodicity. 

In the next section we introduce and describe the 
model, as well as some relevant details of our simulations. 
In section lTTTl we present our numerical results. We report 
on the scaling behaviour of the largest Lyapunov expo- 
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nent with the system size N, and point out connections 
with Hamiltonian dynamics. We analyse the relaxation 
to equilibrium, and we make a systematic characterisa- 
tion as a function of the system parameters. Finally, a 
discussion and summary are presented in section llVl 



II. MODEL 

Our model is a set of N symplectically coupled (hence 
conservative) standard maps, where the coupling is made 
through the coordinates as follows: 



where the superindex T indicates the transposed matrix, 
and J is the Poisson matrix, defined by 



J = 



I 

-I 



(3) 



I being the N x N identity matrix. A consequence of 
Eq. J5J is that the Jacobian determinant \dG/dx\ = 1, 
indicating that the map G is (hyper)volume-preserving. 
In particular, for our model 



dG 
dx 



I I 

B (I + B) 



(4) 



0i(t + l) 
Pi (t + i) 



Oi(t) + Pi (t + 1) (modi), 
Pi (t) + ^sm[2n6 i (t)}+ 



(1) 



— T 



N 



sin[2ir(8j(t)-8j(t))] 



(mod 1), 



where t is the discrete time t — 1,2, . . ., and a > 0. The a 
parameter is the usual nonlinear constant of the individ- 
ual standard map, whereas the b parameter modulates 
the overall strength of the long-range coupling. Both pa- 
rameters contribute to the nonlinearity of the system; it 
becomes integrable when a = b = 0. For simplicity, we 
have studied only the cases a > 0, b > 0, but we expect 
similar results when one or both of these parameters are 
negative. The systematic study of the whole parameter 
space is certainly welcome. Notice that, in order to de- 
scribe a system whose phase space is bounded, we are 
considering, as usual, only the torus (mod 1). Addition- 
ally, the maps are localised in a one dimensional (d = 1) 
regular lattice with periodic boundary conditions. The 
distance is the minimum distance between maps i 
and j, hence it can take values from unity to ^ ( N ~ 1 ) 
for even (odd) number N of maps. Note that r.y is a 
fixed quantity that, modulated with the power a, enters 
Eq. as an effective time-independent coupling con- 
stant. As a consequence, a regulates the range of the 
interaction between maps. The sum is global (i.e., it in- 
cludes every pair of maps), so the limiting cases a = 
and a = oo correspond respectively to infinitely long 
range and nearest neighbours. In our case d=l, thus 
< a < 1 (a > 1) means long-range (short-range) cou- 
pling. Moreover, the coupling term is normalised by the 



sum 



d J x dr r d 



N 1 - 



1-a/d 



-aid 



to 



yield a non-diverging quantity as the system size grows 
(for simplicity, we have replaced here the exact discrete 
sum over integer r by its continuous approximation). 

If G(x) denotes a map system, then G is symplectic 
when its Jacobian dG/dx satisfies the relation 



raipi' 

m 



dx I V dx 



(2) 



where x is the 2iV-dimensional vector x = (p, 6), and 



B 



I K 6l C 2 1 ... C7V1 ^ 
C12 Kg 2 ... C N2 



\ ClJV C2N ■■■ Kg N J 



(5) 



with 



K 9i = acos[27T0i(t)] + — 



and 



b cos[2ir(di{t) -6j(t))] 



N 



where i,j = l,...,N. It can be seen that, 

T 



dG 

dx 



hence 



dG 
dx 



J 



I B 
I (I + B) 



-B I 
-(I + B) I 



(6) 



(7) 



This quantity, multiplied (from the right side) by the 
matrix yields J. Therefore our system is symplec- 
tic. Consequently, the 2A^ Lyapunov exponents Ai = 
Aa/ j A2, A3, X2N are coupled two by two as follows:Ai = 
— A2JV > A2 = — A2Ar_i > ... > An = — Aat+i > 0. In 
other words, as a function of time, an infinitely small 
length typically diverges as e Alt , an infinitely small area 
diverges as e^ Al+A2 ^, an infinitely small volume diverges 



as e 



(Ai+A 2 +A 3 )t 



an infinitely small TV-dimensional hy- 



pervolume diverges as e^*=i Ai ^* (Y]--, Ai being in fact 
equal to the Kolmogorov-Sinai entropy rate, in agree- 
ment with the Pesin identity), an infinitely small (A + l)- 

hypervolume diverges as e^<=i and so on. For ex- 
ample, a (2N — l)-hypervolume diverges as e Alt , and fi- 
nally a 2iV-hypervolumc remains constant, thus recover- 
ing the conservative nature of the system (of course, this 
corresponds to the Liouville theorem in classical Hamil- 
tonian dynamics). 
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For a = 0, similar models exist in the literature though 
in different contexts [l^. l30l l3l| . The present particular 
choice for the coupling was made with the purpose of 
comparison with many-body Hamiltonian systems. In- 
deed, one can derive the map set of equations (1) by 
applying a discretisation procedure to the a-XY model 
with an external field (for more details see GHH). As 
a consequence of having N — 1 terms in the coupling 
summation and the fact that there are N maps, the sim- 
ulation times are of order 0{N 2 ). For this reason it is a 
difficult task to numerically simulate QJ for large values 
of N . To overcome this problem, we used an algorithm 
that takes advantage of the symmetry of the lattice |3^ | 
and shortens the simulation time to 0{N\sxN). 

Initial coordinates and momenta were randomly taken 
from the following uniform distributions: 9^ G [9q — 
69, 9q + 69] and Pi £ \po — 6p, po + dp]. For the coordi- 
nates we used 9q — 69 — 0.5 (i.e., 9i £ [0, 1], homogeneous 
coordinate initial conditions). For the momenta, we con- 
centrated in two cases. To study the sensibility to initial 
conditions we used a uniform distribution over the whole 
phase-space (po = 0.5 and 6p — 0.5). In the analysis 
of the relaxation to equilibrium we used a thin waterbag 
intial condition with p — 0.3 and 6p — 0.05. We also 
checked inhomogeneous initial conditions, in particular 
Po = 9q = 0.3 and 6p = 69 = 0.05 (not shown in this pa- 
per). Note that in this case there is no traslational sym- 
metry in the coordinates. For sufficiently long times, our 
simulations yield the same sensitivity to initial conditions 
that we obtain for the symmetric case 9i G [0, 1]. On the 
other hand, this inhomogeneous initial condition in the 
coordinates has an important influence in the shape of 
the relaxation to equilibrium. This will be further com- 
mented in the next section. We note that the systematic 
study of the role of initial conditions is very instructive, 
but it is out of the scope of the present work. 



III. RESULTS 
A. Sensitivity to initial conditions 

In order to analyse the sensibility to initial conditions, 
we numerically studied the largest Lyapunov exponent 
(LLE) for different values of parameters a, b, a and N. 
We used the well-known method developed by Benettin 
et al [3j| ■ As a consequence of the symplectic structure of 
Q i the Lyapunov spectrum in the 2iV-dimensional phase 
space of the map is, as already discussed, characterised 
by N pairs of Lyapunov coefficients, where each element 
of the pair is the negative of the other. Therefore, the 
LLE sets an upper bound for the absolute value of the 
entire spectrum of exponents. 

We concentrated on the evolution of the LLE for differ- 
ent values of N starting with 6>o = 0.5, 69 = 0.5, po = 0.5 
and 6p = 0.5, Va. In a typical (sufficiently long) re- 
alisation, the finite-time LLE, Am, is a good estimator 
for the analytical definition of the LLE (both quantities 



will coincide when t — > oo 34]). We averaged between 
realisations (typically 100) in order to have small statisti- 
cal fluctuations. We checked that, for appropriately long 
times, \m does not depend on the initial conditions [35j. 



□ a=0.1 




FIG. 1: Lyapunov exponent dependence on system size N in 
log-log plot, showing that Am ~ N~ K( - a \ Initial conditions 
correspond to 9o = 0.5, 66 = 0.5, po = 0.5 and Sp — 0.5. 
Fixed parameters are a = 0.005 and b — 2. We averaged over 
100 realisations. Inset: k vs. a, exhibiting weak chaos in the 
limit N — > oo when < a J$ 1. 

In Fig. n we show the TV-dependence of Am for typical 
values of the interaction-range parameter a, and fixed 
values of parameters a = 0.005, b — 2. Our results show 
that the value of Am vanishes with increasing value of 
N (consequently, so does the rest of the Lyapunov spec- 
trum) as a power-law Am ~ N~ K ( a * > for a j§ 1, and is 
a positive constant for a > 1 (k « 0). In the inset we 
detail function of a. This shows that the map sys- 

tem is weakly chaotic for long-range coupling (Am — > 
when N — > oo), whereas for short-range interactions, Am 
remains positive for all N, meaning strongly chaotic dy- 
namics (as expected [3(|,|37j). Interestingly, this result 
is totally similar to the one numerically measured and 
analytically predicted for the a-XY model [3^ , thus sug- 
gesting equivalent behaviours. Indeed, and as stated by 
Anteneodo and Vallejos, this scal ing is typical of systems 
with couplings of the form 1 jr a 39] . Preliminary simu- 
lations suggest that the fact that the weak chaos region 
extends slightly over a = 1 is an expected consequence 
of finite-size and finite-time effects. 

The dependence of Am with the nonlinear parameter a 
for different ranges of the interaction a is shown in Fig. 
121 We can see that, for a < 1, Am decreases with a 
and saturates for a <C 1. This illustrates the influence of 
this nonlinear term. For increasing a, the sensibility to 
initial conditions raises. On the other hand a has almost 
no effect when a > 1, where the A m is approximately 
constant on the whole a-range. We verified that for a > 1 
a slight increase appears for a > 1 (as in the long-range 
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generically, the sensibility to initial conditions is strongly 
governed by the range of the interactions. Orbits are 
clearly chaotic in the short-range case (as expected) 
and strong numerical evidence of vanishing LLE emerges 
when interactions are long-range, revealing a weakly 
chaotic regime. 

Our present results for the maximal Lyapunov expo- 
nent can be approximatively summarised as follows: 

X M (a,b,N,a)ocf a (b)N- K ^, (8) 

where < a « 1; < b ^ 2; N » 1, a > 0, and f a (b) 
is some function of (6, a) (e.g., f a (b) oc b 63 for a > 1). 

B. Temperature evolution 



FIG. 2: Lyapunov exponent dependence on a for different 
values of a. Fixed constants are N = 1024 and 6 = 2. Initial 
conditions correspond to 6o — 0.5, 89 — 0.5, po = 0.5 and 
Sp = 0.5. We averaged over 100 realisations. 



case), but this effect is negligible. 



As stated above, system is symplectic, hence (hy- 
per) volume preserving, like Hamiltonian systems. A con- 
sequence of this is that 6 may be interpreted as a "co- 
ordinate" variable and p as the conjugate "momentum" . 
We may define a concept analogous to a temperature as 
twice the mean "kinetic energy" per particle |19ll28|. 




FIG. 3: Lyapunov exponent dependence on b in log-log plot. 
Fixed constants are N = 1024 and a — 0.005. Initial condi- 
tions correspond to 8q = 0.5, 89 — 0.5, po = 0.5 and 8p — 0.5. 
We averaged over 100 realisations. Inset: Same data in linear- 
linear plot. 

In Fig. we exhibit the dependence on the coupling 
parameter 6, also varying a. For high values of b (b ^ 1), 
— c(a)6 ' 63 Va, being c(a) a nonlinear function. 
Preliminary simulations show that this exponent varies 
only slightly with a and N, suggesting that it may be 
unique. In the case b — 0, we verified that Am = 0, V<x 
For a > 1 this scaling is valid even for (;« 1. On the 
other hand, if a < 1, a deviation from this power-law 
behaviour emerges when 6^1. 

Our characterisation illustrates the fact that, quite 



rw^E^i-M)) 2 )- (9) 

1=1 

where (...) denotes an ensemble average. This quantity 
can be interpreted as a dynamical analog, and plays a role 
similar to the physical temperature. We refer to as the 
BG-temperature, Teg, the temperature associated with 
an uniform ensemble distribution in phase space. This 
quantity may be analytically calculated, 







/ dPi Pi ~ 




Jo 





which yields T BG = 1/12 ~ 0.083 (VA). 

We studied the evolution of the dynamical tempera- 
ture T for typical values of the parameters as described 
above, focusing on the relaxation towards Tbg- We used 
waterbag initial conditions with homogeneous coordinate 
distribution 9i € [0, 1], and momenta centered at po = 0.3 
with width Sp = 0.05. The manner in which the relax- 
ation takes place depends strongly on the initial condi- 
tions. For example, using < po <^ 1 makes the system 
to temporarily reach a temperature greater than Teg, as 
a consequence of the periodic boundary conditions (torus 
(mod 1)). In general different initial conditions produce 
different relaxation shapes, but scaling with A remains 
similar. We also checked inhomogeneous coordinate dis- 
tributions as initial conditions (8 = 0.3, 89 = 0.05), 
and obtained results qualitatively similar to those corre- 
sponding to the homogeneous case we present here. Al- 
though the temperature value in the metastable state 
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is not as low as in the homogeneous case, long-lasting 
metastable plateaux appear (very much as in the a = 
case of the a-XY model, also known as Hamiltonian mean 
field (HMF) for initial magnetisation M=10), The 
duration of these plateaux also scales with N in a qualita- 
tively similar manner as for the homogeneous case. The 
particular choice of initial conditions used in this work 
yields a rather smooth relaxation with, for example, only 
one inflexion point, thus simplifying our analysis. The 
whole scenario is consistent with the conjecture advanced 
in HJ. 

This type of relaxation has already been reported for 
the particular case a = 0. A two-step process ap- 
pears: firstly a stage where T < Tbg, and then a final 
relaxation to the predicted temperature Tbg- The initial 
regime varies very slowly in time yielding quasistationary 
(QS) states. 
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FIG. 4: Upper panel: Temperature evolution for a = 2 and 
a = 0.6 and four system sizes N = 100, 400, 1000, 4000. Ini- 
tial conditions correspond to Qq = 0.5, 88 = 0.5, po = 0.3 
and Sp = 0.05. Fixed constants are a = 0.05 and 6 = 2. For 
a = 2 the four curves coincide almost completely, all having 
a very fast relaxation to Tbg- For a = 0.6 the same sizes 
are shown, growing in the direction of the arrow. Left bottom 
panel: crossover time t c vs. N, showing a power-law depen- 
dence t c ~ N 0(a) with /3(a) > 0. Right bottom panel: 
vs a shows that for long-range interactions the QS state life- 
time diverges in the thermodynamic limit. Note that when 
a = 0, f3 = 1, and hence t c oc N. Given the non-neglectable 
error bars due to finite size effects, the relation /3 = 1 — a is 
not excluded as possibly being the exact one for < a < 1 
(f3 = otherwise); more precisely, it is not unplausible that 
*c <x ^rr^ (hence, for a = 1, t c oc InjV). 

In Fig. 0| we show the temperature evolution for a — 
0.6 and a = 2 for sizes N = 100,400,1000,4000. The 
curves crossed by the arrow correspond to the a = 0.6 
case, being the first one (from left to right) N = 100. All 
four curves for a — 2 relax approximately at the same 
time, so it has the appearance of a single curve. For a = 
0.6 the typical two-step relaxation is obtained. From now 



on, we define the crossover time t c from the QS state to 
the BG equilibrium state by means of the inflexion point, 
i.e. the time that corresponds to a maximum in the time 
derivative of T. The dependence of t c with N for this 
choice of parameters and initial conditions is reported 
in the bottom left panel. The crossover time scales as 
t c ~ N 13 ^ V a. For a ^ 1, (3(a) f» and then t c remains 
constant (as depicted for a — 2 in Fig. 0J. For a ^ 
1, P(a) > 0, i.e., t c diverges in the thermodynamic limit 
N — > oo. This result indicates the inequivalence, for long- 
range interactions, of the orderings t — > oo and then N — > 
oo, and N — > oo and then t — > oo. Consistently, these 
QS states become permanent (and therefore definitively 
relevant) when N — > oo. Once again, the situation that 
is found coincides with that of the a-XY model 9] . 




FIG. 5: Temperature dependence on a. Fixed constants are 
N — 100 and b — 2. Initial conditions correspond to 6o = 0.5, 
86 — 0.5, po = 0.3 and 8p = 0.05. We averaged over 100 
realisations. 

Finally, we studied the dependence of t c with the non- 
linear parameter a, for different values of a, fixed cou- 
pling parameter b — 2, and size N = 100. In this case, as 
well as in the other situations analised in this work, differ- 
ent behaviours are reported for short- and long- range in- 
teractions. For a > 1, the value of t c , in the limit a — > 0, 
tends to a finite value. The situation is different in the 
long-range case, where, for a <C f , the crossover time t c 
scales as t c = d(a) a 184 being d(a) a nonlinear function. 
Preliminary calculations show that the deviation for the 
case a = 0.9 and a <C 1 is due to finite-size effects. This 
scaling law implies that, for long-range interactions and 
vanishing nonlinear parameter a, the system stays in the 
metastable regime permanently. 

Our present results for the crossover time t c can be 
approximatively summarised as follows: 



t c (a, b, N, a) oc a^^g^b) N^ a \ 



(11) 
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where 0.001 < a « 0.1, < b ^ 2, N » 1, < a < 1, 
and g a (b) is some function of (6, a). 

IV. CONCLUSIONS 

In this work we have analysed a system of N sym- 
plectically and globally coupled standard maps, for both 
short- and long-range interactions. Our results connect 
with studies made in long-range Hamiltonian systems as 
well as with other map systems with vanishing largest 
Lyapunov exponent. We studied a suitably defined dy- 
namical temperature for a region of the parameter space, 
and report the appearance of long-lasting quasistationary 
states, followed by a relaxation to the predicted equilib- 
rium value. The relaxation time scales as t c ~ N^ 1 -"' , 
with (3(a) w when a ^ 1, and (3(a) > when 
< a J$ 1, thus diverging as N — > oo. This is a fea- 
ture also present in the a-XY Hamiltonian model, and 
constitutes a major conjecture in nonextensive statisti- 
cal mechanics (NSM) . Regarding the sensitivity to initial 
conditions, we calculated the maximum Lyapunov expo- 
nent Am as a function of the different system parameters, 
namely, the size of the system N, the nonlinear parame- 
ter a and the coupling constant b. We found that, as the 
number N of maps grows, Am vanishes as N~ K ^ with 
n(a) > for < a < 1, and k w for a ^ 1. The de- 
pendence of k with a is reported and compared to be the 
same as in the a-XY model. We note that, even though 
the initial evolution of Am depends strongly in the initial 
conditions (as the temperature does), this scaling does 
not. In fact, it coincides with that calculated with uni- 



form initial conditions (i.e., initial Pi,9i uniformly dis- 
tributed in the [0,1] interval) for t 3> 1. Our results 
exhibit that, in the presence of long-range interactions, 
the system tends to be weakly chaotic in the thermody- 
namic limit. This is another feature also present in the 
nonextensive theory that suggests its applicability. The 
lack of ergodicity exhibited in our model is thought to 
be related to a (multi)fractal constraint in the available 
phase-space that is enhanced in certain limits such as 
N — ► oo. This possibility would be connected to ex- 
tremely low Arnold diffusion caused by the remanent of 
KAM tori and islands, as has been suggested for simi- 
lar high-dimensional systems \^ . The similarities of the 
present coupled map model and the a-XY Hamiltonian 
model suggest that both share features that may have a 
common dynamical behaviour. This fact, together with 
recent results obtained for low-dimensional maps, may 
clarify the role that NSM plays in the correct description 
of long-range dynamical systems. 
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